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We discuss the nature of the radial and azimuthal components of centrifugal force 
associated with fluid flows in the equatorial plane of black hole space times. The 
equations of motion are solved for the radial and azimuthal components of the 3- 
velocity V which are then used in evaluating the nature of the various components 
of inertial accelerations. It is shown that the reversal of centrifugal force is governed 
mainly by the dominance of the azimuthal velocity and the reversal occurs for r, 
mostly at 2m $ r ^ 3m, depending upon the boundary condition. 
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1 Introduction 

Most of the discussion on fluid flows in astrophysical context are Newtonian S&i in 
nature, with a few considering post Newtonian features □'BEli]. However, if the flows 
are to be considered in the context of black holes, it is only natural to expect the 
effects of space time curvature to affect the flows considered in a general relativistic 
framework. Such considerations may indeed show new features, quite unfamiliar in 
pure Newtonian descriptions. One of the ways to study such features is to develop 
a full scale general relativistic formalism of flows around black holes. However, 
there is one aspect of general relativity which yields the possibility of looking for 
Newtonian features, through a partitioning of space time in a 3+1 language. It has 
been shown by Abramowicz et al. B that one can indeed introduce the concept of 
'inertial forces' in general relativity through a covariant fornialism, as applied to 
test particle motion in curved space time. Recently Prasanna □ (hereafter refers to 
as Paper I), using ajaatching of the 3+1 conforntal splitting of space-time (optical 
reference geometry 113) , with the ADM splitting iA in terms of the lapse function 
a and shift vector /3* has expressed the components of the covariantly defined 
centrifugal and Coriolis forces acting on a fluid element in the given curved space 
time to be 



VV^d, 



V 



and 



where. 



a 
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(4) 



a 



Here, it^ — (u*, u*) is the four velocity of the fluid element and $ is the gravitational 
potential of the back ground space time. Further it was shown in Paper I □ that 
in static space time for the case of pure dust fluid (pressure zero) the centrifugal 
force would show reversal only when is non zero whereas when both and 
are non zero there is no reversal. It is known that generally in most of the 
stationary flows the radial component of the velocity vector is non zero and one 
should see what happens to the force components if the pressure is non zero and 
consequently the gradient of pressure also comes into action. It is possible that the 
gradient of pressure may balance the force due to radial velocity and thereby allow 
the azimuthal component of the velocity to dominate and consequently affect the 
centrifugal force. 

With this in mind, we now consider the study of inertial forces for general 
flows with pressure non zero and with different equations of state. We shall in this 
consider the flows essentially confined to the equatorial plane = -k /2 and evaluate 
the inertial force components, after solving for the velocity fields from the standard 
equation of motion on a given background geometry. 

2 Formalism 

The general equations of motion for a perfect fluid on curved space time are given 



Using the definition of 3- velocity V^' , one can write explicitly the equation of con- 
tinuity 



by: 



(p + P) .y + /t^"P, = 



(5) 




(6) 



) 



.ov dP 



and the momentum equations 




(7) 




) 
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2 



(8) 



We shall presently consider these on a stationary back ground as specified by 
the Kerr geometry in Boyer-Lindquist frame: 

ds^ = --^ ^ -dt^ - 4a ^ dtd<j) + '--dr^ + p^dO^ + ^sin^Odd,^ 9) 

with A = - 2mr + a^, p'^ ^ r"^ + S = (r^ + a^)^ - 

The various components of the shift vectors (3i and the 3-metric jij , alongwith 
the lapse function a are given by (on the plane 9 = 7r/2), 



o + — 2TOr 

A = (0,0,0,-2am/r), 



(10) 
(11) 



2r7i 

7rr = 1 V ^ 



r r- 

while the gravitational potential <& is given as 



and = + + 



2mar 



(12) 



1 



r + a r — 2mr 



(13) 



Using these in the general equations referred to, alongwith the components of 
connection coefficients for the metric (^) one gets, 
the momentum equation: 



,^ .r-r-"- ] 0? — rrir 2m(a? + r'^)^ 2am, , ^ 9x /-r^2 -rA^-, -^ 



(ir 



Ar 



Ar2 



A 



rN2 



+ ir(^^) 



r4 
4am 



-(^^)^ 



Am 
—J 

dP 



du"^ 2am, , „ 9x, ax9 2 , 2ma 



and the equation of continuity: 
1 dvT 2am 

1 



u'' dr l\r'^ 



-^{a^ ^ r'') ^ l-{a^ - mr) 
I\r-^ rA 



Ar^ 



(ma + r (2m — r)) = 



1 



{P + p) dr 



(15) 



(16) 



where, 
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= aV'- I3\ 



(17) 



Wc shall consider the flows with different equation of state and look at the 
behaviour of the force components. For the present we are restricting ourselves to 
the case of flows on the equatorial plane {6 = n/2) of the central gravitating source, 
a condition which has been incorporated in the equations above. 

3 Solutions 

Case (1 ), P^O; Dusty Fluid 

As this refers to the case when streamlines are geodesies, the results are same as 
obtained in the case of test particles, considered in 'Paper F. 

Case (2), p=constant 

(a) Schwarzschild space time 

We shall first consider the static case by putting a = and consequently have the 
equations: 



d 



dr 



2m 



3m 



2m 



1 - 



2m 



1 dP 



r J {P + p) dr ' 



dr r 
and the equation of continuity: 



1 - 



2m 



(19) 



d 



dr ^ r \ r 

2m 



2m 



dP 

{P + p)~di^' 



Thus one finds aV^ = (l — ^) and aV^ satisfies the equation 



(20) 



- 1- 



Solving which one can get 



aV^ dr r 



2m, 
r 



2 5m 



3m 
1 - 



2m 



r J 



8 m, 



0. 



V r / r'^ \ r 



[l + cir4(l- 2M)] 



(22) 
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where ci is the constant of integration to be chosen appropriately. The components 
of centrifugal force are then given by 



P [{2m- r){cir'^{6m'^ - 5rm + r'^) - 1)P + {3m-2r)r'^] 
" r3(r3 + P{2m - r)){ciP{r - 2mf + 1) 

/ [(2rn — r)(cir^(6TO^ — 5rm + r^) — 1)P + (3m — 2r)r-^] 

(C/)^ = , = — ^ • (24) 

"'^''r^'"'"^ - 2m)r3 + l(ci/2(r - 2to)2 + 1) 

It should be noted at the outset that for / = 3^/im, at r = 3m, the total velocity 
V tends to 1, the velocity of light and thus makes the Lorentz factor infinite. This 
feature is independent of the constant ci as may be seen from the expression for V 
(= [l + ciP{r-2mY]/[l + cir'^{r-2m)]). Thus, restricting the value of / < 3^/3m, 
it may be seen that the centrifugal reversal occurs for two values of r one close to 
r = 3m and other close to the horizon r — 2m, depending upon the choice of I and 
ci. For very low values of I and ci the real roots may all lie within the horizon (e.g. 
for I — 1 and ci = 6, ri = 1.73, r2 — 0.175, = —0.153). In principle there are 
five zeros for the quintic equation obtainable by setting Cf = 0, but only two of 
them lie outside the horizon. For a fixed I as ci is increased, the outer most root 
stays fixed at r = 3m, while the other one moves closer and closer to the horizon. 
Figures 1 gives the plots for the components of the centrifugal force for different 
values of I and ci. 
(b) Kerr space time 

The equations of motion for this background with p = constant may be obtained 



from the general equations (|lj), (15) and (|lq). Though the equation for the az- 
imuthal velocity is decoupled from the other two, the equations are not tractable to 
solve analytically. Thus integrating them numerically with the appropriate bound- 
ary conditions one can obtain the profiles for the velocity, using which the profiles 
for the force components may be obtained. Figure 2 gives the plots for the various 
profiles for different values of the Kerr parameter a. As may be seen, there is a def- 
inite reversal of both components of the centrifugal force in each case at locations 
close to r < 3m (see Table 1). The pressure profile shows a monotonic increase 
close to the central compact object. 



Case (3), p = P; Isothermal Fluid 

(a) Schwarzschild space time 

In the extreme case of isothermal equation of state the equations of motion on the 
Schwarzschild background reduce to 



d , TO / 2m \ 3m / 2m 



-1 



-{r -2m){aV'l'f 



2m\ dP 
r I dr ^ 
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and 



ar r 



^{aV^ + - ( 1 
ar r 



2m 
r 



2m\ 



1 _ = 0, 



3™ 



aV =0 



(26) 
(27) 



1 - 



2 TO 



1 



— ) (ay 



rN2 



(28) 



It is seen from (g^) that the angular velocity shows the same behaviour as in Case 
2(a) with aV^ = 73- (l — The radial velocity comes out constrained from 

the continuity equation equation to be aV^ = 75- (l — ''^i being constant of 
integration, which by dimensional considerations and also due to the consistency 
with the equation of continuity can be identified with the ratio M/p, M being the 
accretion rate and p the density. With these, the components of centrifugal force 
turn out to be 
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Fig. 2 : p=constant, Kerr case. It shows the variation of (a) , (b) u'^, (c) Cfr, (d) 
Cftf, as functions of radial distance in unit of mass of the black hole {r/m). Solid, dotted 
and dashed curves indicate the results for Kerr parameter a = 1, 0.9 and 0.5 respectively, 
(e) Variation of iT, u'^, Cfr and Cf^ as a function of r/m for a = 0.9 which shows the 
relation of the centrifugal reversal location with the location where u'^ crosses over u''. (f) 
Variation of pressure P as a function of radial distance r/m. 



{Cf)r = 



7^ 



2m 
r 



2vl 

„2 



3m 
r 



(29) 



{Cf)4 



J.2 I ^ 



(30) 



with 7^ = (1 - V"^) ^ and V"^ = V^Vr + V^V^. From these expressions it is again 
clear that these do not change sign for any r > 3m. However, for r < 3m, both 



components of the centrifugal force change sign for r 



3 m 
2 



1± \/l 
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Fig. 3 : P = p, Schwarzschild case. It shows the variation of (a) , (b) V^, (c) Cfr, 
(d) Cf^ as functions of radial distance in unit of mass of the black hole {r/m). Solid, 
dotted, dashed and long-dashed curves indicate the results for l/v\ = 2, 1.5, 1 and 0.9 
respectively, (e) Variation of V, , Cfr and Cf^, as functions of r/ra for l/vi = 2, 
which shows the dominance of azimuthal velocity over radial velocity and the reversal of 
centrifugal forces, (f) Variation of pressure P as a function of r/m. 



Obviously the reversal depends upon the reality of r, which demands — < 
Figure 3 shows the profiles of velocity and centrifugal forces for four different values 
of the ratio — . For ^ < 1 the root lies inside the horizon r = 2m and then the 
reversal can not be seen. On the other hand for — > 1 the root lies outside the 
horizon r = 2m but at the location < 3m. 

(b) Kerr space time 

In the case of Kerr background, again the equations are integrated numerically 
and the profiles presented in Fig. 4. The plots give again the velocity and force 
components alongwith the pressure profile for different values of a. Here again as 
in the case of p constant, the centrifugal force changes sign at a radius r between 
2m and 3m, the exact location of the zero depending upon the boundary condition 
(see Table 2). 
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Case (4), P = ; F being the gas constant 

The adiabatic equation of state yields for the static space time, the fluid equations 
which are coupled for V"^ and P, while the azimuthal component of the velocity 

again has the same profile = ^ (l — "^Y^^ ■ The equations for V and P are 
solved numerically for diS^erent values of F and the results are presented in Fig. 
5 for three different values of F; 5/3, 4/3 and 1.0001. The radial velocity profile 
shows difi^erent slopes for different values of the gas constant 

which behaviour is reflected in the profiles for pressure and density too. The az- 
imuthal velocity profile is independent of F. Again the centrifugal force shows 
reversal in the region 2m < r < 3to, in accordance with other cases. Though 
both components of the velocity show a decreasing trend, calculations show that 
the total velocity V tends toward c(= 1) as r ^ 2m, which is in accordance with 
expectations. 

Figures 6 and 7 show the profiles for the same physical quantities in the case of 
Kerr background for the two distinct eases, F = 4/3 and F = 5/3 and for different 
values of the Kerr parameters 'a'. The behaviour of various quantities are just 
similar, again showing the feature of centrifugal reversal in the region 2m < r < 3m 
which is outside the ergo region. 
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Fig. 5 : Schwarzschild case for adiabatic equation of state. It shows the variation of (a) 
, (b) V^, (c) C/r, (d) C/0, (f) P as functions of radial distance in unit of mass of the 
black hole (r/m). Solid, dotted and dashed curves indicate the results for gas constant 
r = 5/3, 4/3, 1.0001 respectively, (e) Variation of V , V'*' , Cfr and CU as functions of 
r/m for F = 4/3 which indicates the dominance of azimuthal velocity over radial velocity 
and the reversal of centrifugal forces. 



4 Discussion 

We shall first consider the azimuthal component of the centrifugal force in the 
Schwarzschild geometry for the case of dusty fluid, which turns out to be 



(Cf) 







I 

2^ 



2mc 



1 



2m\ 



1/2 



(31) 



|and shows no reversal just as the radial component as considered earlier in Paper I 
□. On the other hand when we consider the Kerr background the force components 
do show a reversal even when P = 0. Obviously the extra contribution towards 
angular velocity comes from the dragging of inertial frames, an effect that dominates 
as r ^ 2m the ergo surface. Thus one finds the reversal of the centrifugal force to 
occur close to the location where u'^ crosses v/". 

When the density is constant with pressure non zero, the analytical expressions 
are obtained for the components of the centrifugal force in the static case, which 
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Fig. 6 : Kerr case for adiabatic equation of state with F = 4/3. Plots are similar to the 
ones as Fig. 3. 



show that though in principle there can be five roots for the equation (C/) = 0, 

only two of them are outside the event horizon r = 2m. However the realisation of 
these two roots also depend upon the condition that I be < 3-\/3to. For I = SV^m, 
V ^ c, as r ^ 3m, which indeed is interesting, for the minimum impact parameter 
for a zero rest mass particle in the Schwarzschild space time is B-v/Sm.For this value 
of I the centrifugal force is infinite at r = 3m since the Lorentz factor goes to 
infinity. However, for other values of I the components have two zeros, one close 
to 3m and the other close to the horizon. The occurrence of the two zeros in the 
region outside the horizon does indeed depend upon the choice of the constants / 
and ci. The reversal for a given I, if occur at r = 3m remains unchanged for all ci 
whereas the one nearer the horizon gets closer to the horizon with increasing values 
of ci . It is also important to note that for no value of I and ci the reversal occurs 
for r > 3m. 

In the case of Kerr space time, where the equations are integrated numerically 
and one needs to choose the boundary conditions appropriately, it may be seen that 
the location of the reversal does change with the relative magnitudes of the velocity 
components at infinity. Table 1 gives the locations of the reversal for different values 
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Fig. 7 : Kerr case for adiabatic equation of state with F — 5/3. Plots are similar to the 
ones as Fig. 3. 



of 'a' with corresponding boundary conditions for p=constant case. As may be seen 
in the first case, > (the subscript 'o' refers to the values at boundary from 
where the integration starts), and this is maintained throughout the flow with the 
azimuthal component dominating, the centrifugal force reversal occurs at the same 
location as for a test particle Bll3. On the other hand in the second case though 
> uf, as the flow approaches the black hole the u'^ takes over u'^ at around 
r = 20m leading to the reversal of centrifugal force at the location for example 
r ~ 2.783m for a = 0.9 (Fig. 2e), a location slightly closer to the black hole than 
the reversal radius of 2.822m for the test particle case. The other root seems to 
occur well inside the horizon and thus of no consequence. 
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Table 1 
p=constant 



a 





0.5 


0.9 


1 


Boundary Values 












ul = 10-« 


R 


3.0 


2.944 


2.822 


2.783 


ut = 10-6 












Po = 10-5 


Ri 


2.955 


2.904 


2.783 


2.744 


= 5 X 10-5 


R2 


2.117 


1.948 


1.596 


1.371 


ut = 10-6 












Po = 10-3 



When we consider the isothermal equation of state (P = p) in the Schwarzschild 
background one could solve the equations analytically yielding the velocity compo- 
nents in terms of two constants vx and I. As mentioned earlier v\ can be identified 
as the ratio M/p and it is seen that the expressions for the components of cen- 
trifugal force depend upon both vi and /, with the reversal occurring at a value of 

r = ^ ^1 ± \/ 1 — 9F^^ • ^^'^ "-^^ straight away see from this expression that 

the reversal would occur at r = 3rn for < 1, which is in accordance with earlier 

observation as in the case where would dominate over Tf . On the other hand 
when ^ > or equivalently > then the centrifugal will never be zero for 
any r and thus no reversal would occur. This again confirms the fact that when 
the radial velocity dominates over the angular velocity there is no reversal of the 
centrifugal force. Table 2 presents the location of centrifugal reversal for the case 
P = p, which distinctly shows the possibility of the reversal occurring either at the 
same position as for a test particle or at locations different but inwards and closer 
to the black hole. 



Table 2 



a 





0.5 


0.9 


1 


Boundary Values 


R 


3.0 


2.944 


2.822 


2.783 


ul = 10-» 
ut = 10-6 
Po = 10-5 


Ri 
R2 


2.992 
No root 


2.936 
No root 


2.813 
1.438 


2.773 
1.058 


uj; = 5 X 10-5 
ut = 10-6 

Po = 10-3 



In the case of Kerr black hole one set of locations are always outside the ergosurface 
r = 2m, while the other could be inside. 

Coming to the case of polytropic equation of state Table 3 presents the locations 
of the centrifugal reversal for different values of a and of gas constant T (= 1 -l- ^, n 
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being the polytropic index). It is interesting to see that for F > 5/3 the location of 
centrifugal reversal coincides with those of a test particle, indicating the dominance 
of u'^ over u^. On the other hand for T, 4/3 < F < 5/3, the reversal occurs at 
locations slightly higher than those corresponding to the test particle case. As may 
be seen from Fig. 6e and Fig. 7e, the location of the intersection of u"^ and u*^ 
profiles 



Table 3 



F 


a 





0.5 


0.9 


1 


n 


4/3 


Ri 


3.068 


3.017 


2.90 


2.860 


3 




2.005 


1.872 


1.447 


1.091 




1.45 


Ri 


3.011 


2.958 


2.838 


2.799 


- 2.22 






2.003 


1.870 


1.444 


1.081 




3/2 


Ri 


3.005 


2.950 


2.830 


2.790 


2 






2.003 


1.870 


1.443 


1.078 




5/3 


Ri 


3.0 


2.945 


2.823 


2.783 


1.5 




R2 


2.002 


1.869 


1.442 


1.072 




2 


Ri 


3.0 


2.945 


2.823 


2.783 


1 




i?2 


2.002 


1.868 


1.441 


1.067 





seem to have a direct bearing on the location of the centrifugal reversal. Figure 8 
shows the variation of the location of centrifugal reversal as a function of the gas 
parameter F, for different values of the Kerr parameter a. 

If we consider the Coriolis force as given by (H), it is seen immediately that 
while for static space time it is zero, as it should be, in Kerr space time only the 
radial component is non zero in Boyer-Lindquist coordinate system and is given by 



ar + o?r + 2ma'^ J 

which matches exactly with the expression for the test particle case. This in fact 
is to be expected as Coriolis force does in no way depend upon the radial velocity, 
in the present context. 

In conclusion we can say that for a fluid in curved space time, the centrifugal 
force reversal occurs only through the dominance of azimuthal velocity over the 
radial velocity. For Kerr black hole, due to the effect of inertial frame dragging, 
which always goes to increase the effective angular velocity, the reversal can occur 
even when the pressure of the fluid is zero. On the other hand when pressure is 
non zero, the gradient of pressure seems to bring down the effect of radial velocity 
and allow the fluid to have larger azimuthal velocity and thus lead to centrifugal 
reversal both in radial and azimuthal directions. The above discussion, however 
did not take into account any viscosity, which if present may affect the azimuthal 
velocity through the transfer of angular momentum and thus increase the effect 
of radial velocity. It would be necessary to extend these studies by including the 
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Fig. 8 : Variation of the location of centrifugal reversal as a function of F for Kerr 
parameter (a) a = 1 (solid curve), (b) a — 0.9 (dotted curve), (c) a — 0.5 (dashed curve) 
and (d) a = (long-dashed curve) . 



viscous forces and analyse the effect on centrifugal reversal. 
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